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sutuw 

Scat  an»>diJMn«lauLl  fission  proossssa  art  eonsldtrsd  vlth  par- 
tleulsr  rsfsrd  to  qusstloas  of  erltlesl  lsz>cth  snd  SDsrgy  dspsodence 
of  tbs  products  of  flaeioe.  Various  lots r conns ctlosis  betvssn  the 
functions  ussd  in  tbs  invariant  ibbsddlnf  spproaeb  aikd  tbs  classical 
approaob  arc  Indleatsd.  In  particular  it  is  sbovn  bov  one  asy  pass 
frcB  tbs  eonsidsrstion  of  csrtaln  linear  tvo-polnt  boundary  valus  problsas 
to  Donlinsar  (KlceatiOB  type)  initial  valus  problssw.  lbs  forasr  arc 
ussful  for  tbsorstloal  eons  ids  rat  ions,  tbs  lattsr  for  ntasrleal  eas^utation. 


D2FIiri?I0S3  OF  VAaiABLgS 


r-a2^2 

il 


9(n»x)  •  tbs  probeblUty  that  n  aautroos  ar«  reflaetad  over  all  tlae 
by  a  rod  of  length  x  aa  a  result  of  aa  loeldent  neutroo. 

“<*»**)  -  **" 

0(*)  -  np(ii.x) 

DaO 

-  the  probahlUtj  that  a  particle  starting  to  tbs  left  at  j  at 
tlae  saro  viu  result  In  n  partleles  erglng  fro*  over  all 
tlM. 

qj,(B,7)  -  the  seas  probablUtsr  asstalng  that  the  particle  starts  to  the 
right  at  tlae  zero. 

s 

fir, 7)  - 

n-O  ^ 

■  g  9R(»#y)  r“ 

My)  -  £  o^T  (n,y) 

ttol  *• 

V(x)  -  the  expected  mM>er  of  neutroas  eaergent  fna  oae  end  of  a 
rod  of  length  x  as  a  result  of  a  trigger  neutron  Ineldent  at 
the  other  sad. 

e 

''g^***)  ■  expected  sadwr  of  neutrons  over  all  tins  passing  t  toeard 
the  right  as  a  zusult  of  one  tz’lgger  neutron  entering 
bar  [x,0  ]  at  x. 
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v^(b;z)  -  tha  •orraspoodlng  mabar  paasloK  to  tba  laft. 

*  tba  axpaetad  aiotwr  of  nautrona  ovar  all  tlaa  aovloc  to 
tha  rlsbt  at  tha  point  t  aa  a  rasult  of  a  aourea  nautroc 
■ovlag  to  tba  rlgbt  at  7,  tba  lanctb  of  tba  bar  balng  z. 

V(z,t)  ■  tba  azpaetad  mabar  of  partlelaa  aarglnc  froa  x  batvaan 
tlaa  zero  and  tlaa  t  froa  a  rod  of  laacth  x  dua  to  a 
tri^tar  navtron  antarlng  at  z  at  tlaa  zaro. 

u(z,t)  -  U^(x,t) 

T(x,t)  -  tba  eorraapoadli^  mater  traaaitted. 

■  tha  azpaatad  aaiber  of  nautrona  aaarglD4  froa  x  ovar  tha 
tlaa  Intarval  [  0,t]  dua  to  an  Initial  naxitron  aovlng  to 
tha  left  at  7  at  tlaa  zero. 

Up(7*t)  *  tha  eorraapoadlng  nuAar  due  to  a  neutron  aovlng  to  tha 
right. 

u(z,v)r)dr  -  tba  ezpaetad  aihber  of  nautrona  rafleetad  froa  0  vlth 
energ7  batvean  r  and  r  >  dr,  over  all  tlaa,  dua  to  an 
Initial  nautron  of  energj  v  entering  the  rod  at  0  at  tlaa 
zero. 

T(z,v;r)dr  >  tha  expected  madwr  tranaaltted. 

q(a,r)dr  -  the  probahUltp  of  energp  ebai^  froa  level  a  to  a  level 
batvean  r  and  r  dr  due  to  aflaaionlng. 

«l|(7«v;r)dr«  tba  ezpaetad  m^wr  of  nautrona  aaBigliig  froa  0  vlth 
eaarg7  batvean  r  and  r  -f  dr,  over  all  tlaa,  dat  to  an 
Initial  neutron  at  anarg7  w  starting  to  tha  right  at 
7  at  tlae  zero. 

eorreopondlng  qunntit7  when  the  neutron  atarta  to 


the  laft 
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ruBcnoiuL  BQUAnons 

Richard  Ballaas  and  Robert  Kalaba 
nie  BAUD  Corporation,  Santa  MOBlea,  California 

and 

0.  Nlltoo  Win« 

Loe  Alaaoo  Sdeotlfle  Laboratory,  Uhl varsity  of  California 
Los  Alans,  Itav  Italeo* 

I.  Introduction 

In  our  prevlou5  papers  on  the  applleatlOD  of  the  principle 
of  invariant  lid>eddin£  to  neutron  transport  prooesses,  [  3]  ,  [  ^  ]  ,  se 
considered  the  prbblosi  of  deterainlne  probability  dlstrlbutloDS  for 
the  reflected  and  transBltted  fluxes  associated  vlth  one-dlasnslonal 
irods.  Using  only  expected  values,  these  techniques  led  to  asthods  for 
detemining  critical  length. 

'Bae  custoaary  aethod  at  treating  these  questions,  based  upon  the 
equation  of  transfer,  leads  to  linear  differential  equation?  for  the 
expected  values,  vlth  tvo-polat  boundary  eonditicos.  The  erltleal 
lengtl  i  s^aaculated  as  the  solution  of  an  eigenvalue  prbblan. 

The  aethod  presented  in  ^  s]  aiad  [ 4 ]  is  of  an  entirely  different 
nature.  The  basic  equations  are  now  nonlinear  differential  equations 
vhoee  solutions  are  detendned  solely  by  initial  values.  The  critical 
length  turns  out  to  be  the  first  value  of  which  the  solutions  of  the 
equations  for  expected  flux  becoas  infinite. 

- « - 

Work  perfoxBsd  in  part  lasdcr  the  ausploes  of  ths  O.S.  Atoaie 
totrgy  Coadasion. 
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Xb  Tltfw  oC  ttm  fact  that  tbccc  two  a^roatikoa  to  prohlaaa  of  thla 
type  apvMT  to  ho  so  dioparato,  it  Is  of  sobs  latorsst  to  Ilok  thsa, 
and  to  adsrstaad  ths  rslatlOD  hstimoB  tbs  various  functions  that  art 
iatrodassd  ia  thsss  fashions.  In  [  l]  ,  whsrs  Invarlanes  princlplss 
aro  iMod  to  treat  the  closely  related  prooesees  of  radiative  transfer, 
rnlisniBlnn  states  that  V.  V.  Sobolev  has  obtained  various  econeetlocis 
betwaea  ths  solutions  of  ths  equations  derived  froa  Invarlanoe  prinoiples 
sad  the  fwstions  arising  frea  the  classical  equations  of  aathsaatioal 
physios.  Vs  have  been  aahle  to  secure  these  papers,  and  although  we 
have  searehed  the  literature  for  papers  on  this  tople,  we  have  not  been 
able  to  find  any. 

It  is  this  interrelation  which  ve  propose  to  discuss  in  this  paper 
in  SOBS  detail.  We  shall  show  that  the  equation  of  transfer,  as  well 
as  the  fiastional  equations  obtained  prsviously,  aty  all  be  derived 
fr«B  a  MHBon  point  of  view.  In  sobs  eases,  this  Bsthod  yields  linear 
eqwticBS  with  two-point  boundary  eondltions;  in  other  eases,  it  yields 
Bonlinsar  eqwticns  with  initial  value  condltlans;  ia  still  otbsr  eases, 
we  bbtaln  the  desired  relations  oenneetlng  the  solutions  of  the  two 
types  of  equations. 

tSB  great  analytic  utility  of  linear  equations  for  a  rigorous 
study  of  the  character  and  nature  of  diffusion  processes  nuat  oertainly 
be  BVhasised,  ef .  Taller,  [6  ]  ,  Lehner  and  Wing,  [9]  •  Honetheless , 
in  eonBBctlen  with  the  eo^potatlonal  detexmlnation  at  the  solution, 
two-point  bOBSdary  eendltlons  cause  serious  difficulties,  ef.,  for 
niseis.  Davison,  [5  ]  >  ^or  a  dlscuseion  of  a  mmbar  at  Bsthods 

thsae  petblsBS.  Ihe  nonllnesr  equations,  although 


utilised  to  treat 
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wrtainljr  not  &•  handy  for  analytlo  purpoooa,  po«M«a  th*  aarlt  of 
rcqolrloc  only  initial  eonditiona  for  thalr  nxatrleal  solution.  'I^ioy 
sees  veil  suited,  In  eonsequeaee,  for  aodem  dlcltad  eoq^utem.  For 
refereaoes  to  related  work  of  Rsdbeffer  and  Luneherg,  see  [2]  .  In 
eonaeetlon  vltfa  the  traamlsslon  line  equations,  the  paper  of  Fierce, 
[lO"]  ,  Is  relevant. 

Aiqr  rlgox'oos  discussion  of  the  stochastic  processes  ve  consider 
here  would  take  us  too  far  afield.  We  refer  the  reader  Interested  In 
these  natters  to  the  forthconlag  nonograph  by  T.  S.  Harris,  [8]  ,  and 
to  his  expository  paper,  [7]  .  Slnllarly,  to  avoid  digressions  which 
would  divert  us  froi  our  principal  object,  that  of  Illustrating  the 
appUeabllity  of  a  variety  of  functional  equation  techniques,  we  shall. 
In  the  aaln,  consider  only  expected  fluxes.  To  detemlne  full 
probability  distributions  would  require  for  the  general  case  a  study 
of  characteristic  functionals,  a  topic  of  s<ae  coaplezlty.  Only  In 
the  tlne-and-energy- Independent  ease,  do  we  discuss  the  collets 
distribution. 

We  shall  consider  only  one-dinsnslonal  processes,  treating  In 
turn  the  tine- Independent,  energy-independent  ease;  the  tine-dependent, 
energy- Independent  ease;  and  the  tine-independent,  energy-dependent 
case.  In  the  study  of  the  tine-dependent  case,  we  are  gratefxil  to 
T.  B.  Harris  Ibr  the  suggestion  to  Introduce  nosnnt  density  functions 
and  for  nueh  other  helpful  erltlelsn. 


See  M.  Bartlett,  ^  Introductl«  to  Stochastic  Processes ,  pp.  78-83, 
where  aany  references  to  t.<«  wtwA  e/niaVinMie  eay  be  found. 
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1.  A  ifctb— tieal  Nod»l  of  0— -Pl—Mlopal  »«iitrop  Trminport. 

L*t  os  aov  daflM  vhat  w«  shall  assn  bjr  &  oos-dlasnsloiial  nsxttron 
transport  proosss.  A  nautron  Is  eonsldcrsd  to  b«  a  point  aovlac  on  a 
Has  Intsnral  «SDtil  an  svost  occurs  which  causes  It  to  transfora  Into 
two  neutrons  of  slnllar  nature,  ons  nowlng  to  the  left  and  one  to  the 
rl^t.  this  event  will  be  called  fission  and  the  proeess  called  a 
fission  process. 

The  probability  that  a  neutron  vill  undergo  fission  In  an 
Inflnitesinal  Interval  of  length  y  is  taken  to  be  X  y  ♦  o(y),  and 
the  probability  that  it  can  traverse  an  interval  y  vlthout  this 
occurring  is  1  -  ^  y  ♦  o(y) . 

Here  X  represents  a  ’nean  free  path. '  We  have  adjusted  our 
notation  to  confoxa  to  that  current  In  physical  circles. 

fhs  prooess  can  be  oonplieated  by  the  ass»vtlon  of  Inhoaogsnelty, 

1. e.,  X  y  replaced  by  the  fwotlon  1/x  (y);  by  the  Introduction 

of  a  distribution  fwietion  for  the  m^ber  and  direction  of  the  particles 
produced  by  fission;  and  by  taking  aecoiat  of  absorption  and  eoUlolons 
which  asrely  reverse  direction  without  Inducing  fission.  Slnoe  these 
aore  rsalistle  features  asrely  add  arlthastle  coaplexlty  vlthout 
analytic  insight,  we  have  Ignored  their  existence  In  this  prellalnBry 
trsstasnt . 

We  shall,  however,  discuas  the  case  In  which  the  neutron  and  Its 
progeny  posess  energies  which  affect,  and  are  affected  by,  fission. 

2.  Invariance  Principles. 

As  we  have  explained  la  earlier  papers ,  cf.  [s]  ,  our  ala 

is  to  introduce  certain  parwasters  with  ths  pi^opei'ty  that  the  neutron 
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transport  process  described  above  can  be  considered  as  a  sequenoe  of 
stochastic  processes  vhlcb  Induce  transfornatlons  of  these  paraaeters. 
Purthexvore,  It  Is  desirable  to  have  these  paraiwters  represent 
physleally  ■eanlagful  quantities  such  as  tlae,  energy,  length,  and  so 
on,  essentially  qnantltlee  which  correspond  to  ■easurablee  and 
c^aervablas.  This  is  aerely  desirable,  but  should  not  be  Insisted 
t^on  to  the  detrlnent  of  the  analysis.  That  there  are  no  indlspenaable 
variables  la  quite  evident  by  now. 

In  applying  Invarlanee  prlneiplea,  we  can  bonw  sane  terelnolOGr 
free  ooHMSleatlon  theory  consider  each  physical  process  to  consist 
of  three  principal  eleeenta,  a  source,  an  observer,  and  a  tran^lttlng 

or  ' co^Munlcatlon  channel.'  As  we  shal 1  see  below,  the  classical 
eqtmtlons  of  ■athaeatleal  physics,  of  which  the  equatice  of  transfer 
1«  typical ,  arise  when  f tmetlcnal  .equation  techniques  are  applied  to 
processes  defli»ed  in  teii»  of  variables  dascrlblns  the  observer,  or  the 
source.  On  the  other  band,  the  new  equations  of  Invariant  Ihbedddlng 
arise  when  the  process  Is  described  peurtly  In  teiaa  of  the  nedliai  and 
partly  In  terns  of  the  observer  or  source. 

By  describing  a  single  process  partly  In  one  set  of  terns  and 
partly  In  other  teimlnology,  we  shall  obtain  the  desired  equations 
linking  the  various  fomulatlons.  There  Is  nueh  to  be  done  In  this 
direction,  with  this  paper  nerely  a  single  step  in  the  direction  of 
a  ncre  cavrebensive  tlnory.  We  sheOl  restrict  ourselves  hers  to 
one-dlnenslosial  processes. 
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II.  Tlt-Ia<«p«nd«»t.  gaTKy-Indw>flnt  Ca— 

3.  yir«t  FoCTulatioa. 

Ve  first  rsesU  mm  rMolts  darlrsd  la  [  3  ]  . 

OoMldir  a  flalta,  tkOMcnsoua  rod,  ttao  lotomd  frta  z  to  0  as 
ladleatod  bolow. 


-I 

0 


Figure  1 

Lot  tM  start  with  the  ease  where  a  'trigger'  neutroo  enters  toe 
rod  froB  the  left,  at  x,  at  tlae  Mro.  We  shall  aasMs  that  energy  Is 
■sgleeted,  and  rid  oorselves  of  tlM  eonsldoratloae  by  asking  for  the 
respeetlTS  prohehiUtlos  that  prsdseljr  n  neutrons  are  reflected  froa 
the  rod  (vhleh  Is  to  say,  snigi  froa  x  to  the  left),  over  all  tlM, 
for  n  *  0,  1,  2,  *•*  * 

Ihls  prohahlUty  will  be  denoted  by  p(n,x).  In  latroduelsg  this 
notation.  In  plane  of  uhst  euetoMrlly  night  be  called  p(a),  ve  am 
telibemtely  ««»THae  attention  to  ths  fact  that  the  length  of  the 
rod  is  to  be  treated  In  this  section  as  n  essential  pamnster  of  the 
process. 

We  am  thus  l^>oddlng  this  particular  process  In  the  fanlly  of 
prooesoes  of  slallar  natum  in  vhleh  the  rod  can  assiM  any  positive 
length. 

Using  Invert anee  prlnolples,  as  described  ln[.kj  in  detail,  ve 
obtain  the  following  bjeten  of  nonlinear  dlffemntlal  equations. 
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p'(0,x)  »  -  p(O^) 

(3.1) 

p’(»,x)  -  -  (n+l)  x'SCn.x)  ♦  x"S(b-1,x)  ♦  p(n-k,x), 

k-1 

▼alld  for  X  2  0,  aloog  wltb  the  initial  ean&ltlan  p(a,0)  > 

Die  study  of  these  equs^lacs  Is  as  usual  greatly  facilitated  by  the 
introduction  of  the  gsoeratlng  fwurtlon 

(3.2)  h(x,r)  «  ^  p(n,*)  r“  . 

Using  the  eqxiatlons  In  (5.1),  ve  readily  see  that  h(x,r)  satisfies 
the  q\iasl-linear  partial  differential  equations 

(3.3)  h^  =  X*^  (r-l)  h  ♦  X*^  r(h-l)  . 

Either  frta  this  equation,  or  directly  froa  (5.1),  (also  see 
Section  7),  ve  derive  the  Siecati  eqxiatlon 

(3.4)  U'  o  (l40^) 

for  the  expected  maber  of  reflected  neutrons,  U(x).  This  equation  is 
valid  for  0  <  x  vhere  x^  is  the  crltleA)  length  irtiieh  turns  oxit 

to  be  X  . 

As  we  shall  see  belcw,  this  equation  can  be  derived  froa  first 
principles,  without  the  introduction  of  higher  aoasnts,  or  use  of  the 
generating  function.  On  the  other  hand,  it  is  alao  possible  to  obtain 
the  partial  differential  equation  for  the  generating  function  by  these 
direst  teshnlqaes. 
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a»OOOA  FOfBlKtiOP 

la  this  foiaulstioe  at  s  flssloo  proosss,  ws  bsgln  irlth  s  partiels 
sltwtad  «t  s  poiat  y,  latsrloar  to  ttas  rod,  at  tlas  soro  and  starting 
off  sitter  to  tte  Isft  or  to  tte  rl«bt. 

-♦ - ( 

7  0 


K 

X 


ngsrs  2 

Lst  us  dsTlns  tte  protebilltlss 

gj^(n,y)  -  tte  proiteMllty  that  a  partiels  starting  to  tte 

left  at  y  at  tlas  tsro  vlU  rssult  la  n  partlelss 
sasrglag  frcai  x  ovsr  all  tlas. 

(‘‘.1)  <jjj(n^)  -  tte  SMS  protehUlty  assailag  that  tte  partlsls 
starte  to  tte  right  at  tte  Initial  Instant. 

Obssrvs  that  vs  havs  now  s\vprssssd  tte  x-dspandsnes,  Indl  eating 
that  vs  vlsh  to  teop  x  fizsd  throughout  tte  snsulng  analysis,  and  havs 
foeussd  our  sttsstloa  i«cd  y.  Ws  vlsh  to  Uribsd  this  proosss  vlthin 
tte  fMlly  of  prossssss  at  this  type  in  vhiai  y  takss  on  any  valus 
bstwssa  0  aad  x. 

TO  obtain  sqmtloas  eharastsrlzlng  tte  fxmetions  vs  haw  Just 
iBtrodaosd,  vs  follow  tte  path  at  tte  partiels  fna  y  to  y  - 
dopsadsnt  1900  tte  Initial  dlrsetloo.  If,  for  sxaivle,  tte  parti ele 
suffsrs  a  eolllslon  In  [  y,y  -  A  ]  ,  vs  obtain  a  nsv  prsosss  In  vhleh 
eas  partiels  starts  fros  tte  right  at  y  -  &  and  tte  otter  starts  f rca 
tte  left  at  y  -  A;  if  tters  Is  no  soUlslon,  vs  ara 


Isft  vltb  a 
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proo*«8  of  tlM  MW  tyv*  far  an  Initial  partlela  at  7  -  A  rtartln« 
to  tfaa  rlgbt. 

TStking  aeeouDt  of  tteM  posalbllltlas,  it  Is  saay  to  saa  that  va 
obtain  ttaa  ralatlons 


qj^(0,y)  =  0 

(1^.2) 

if® 

qj^fey)  “  qR(n-k,y)  i 

♦  (1  -  a)  qj^(n,y+A)  ♦  o(a), 

for  0<  y*  X,  and  n  -  1,  2,  ...»  with  the  boundary  conditions 
(^.3)  qj^(n»r)  =  /jjj  • 

We  obtain  a  ■iailar  set  of  recursion  relations  for  the  sequence 
^qj^(n,y)^  with  the. boundary  conditions 

qj,(a»0)  “‘^On  • 

yroB  these  appraxlnate  relations,  ve  derive  the  following  nonlinear 
dlffereotial  equations 


q£(n,y)  ^q^t-^'y)  qp(Q-Ay), 

(‘‘.5) 

q^(n>y)  -  qp(n,y)  +  X  ^  5ZqpU»y) 


for  0<  y  <  X. 
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Bm  gOBratlac  fWMtlana,  dtflJisd  by  the  MrlM 

f(T,j)  -  £  q.(B,y)  r“, 

B-O  ^ 

•(r,y)  -  <iA^pT)  r“, 

B-0  ^ 

Mtlafy  th«  nonJlTMiir  partial  diffaraittial  •<iuatlons  aad  boundary 
eoadltlOBS 

-  f(r,y)  f(r,y)  ♦  x"^  f(r,y), 

(^•7)  gy  -  f(r,y)  f(r,y)  «(r,y), 

f(r,x)  -  r  *(r,0)  -  1. 

What  1*  aaaalag,  and  a  graat  advantaca  of  this  forBulatlon,  la 

that  tha  MMta  %(n,r)  and  £  <i^(c,y),  k  -  1,  2, 

i>0.  ^  n-l  ^ 

aatlafy  linear  differential  equatione.  Introduelng  the  expected  values 


u,(y) 


B<lp(n,y), 


uj^(y) 


£  nq^(n»y), 

n-1  ^ 


ve  obtain,  throuch  differentiation  of  eq.  (1».7),  the  linear  systea  of 
differential  equations 


v>p(y)  -  «i.(y)» 

»£(y)  - 


Vy),  0<  y  <x. 
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Tba  booDd&ry  eondltiaiis  art  aow  tvo-point: 

Ujj(O)  -  0, 

('*.10) 

u^Cx)  .  1. 

Critical  Maaa. 

Od*  of  the  principal  alaa  of  a  theory  of  noutroo  tranaport  la  the 
drtaimlnatlon  of  critical  aaaa.  We  equate  the  phyaioal  phenoawnon  of 
criticality  vlth  the  aatbeaatleal  concept  of  alngularlty.  In  both 
analytic  approaches,  ve  seek  a  length  of  rod  vhi^  reaulta  In  an 
Infinite  expected  value,  cf.  [s]  and  ]  . 

Tbe  details  of  the  analysis,  however,  are  quite  different  In  the 

two  f oraulations . 

In  the  first  approach,  we  have  an  ordinary  differential  equation 
of  Rieeatl  type,  with  a  given  initial  value,  which  we  Integrate 
mJBsrically  point— by— point  until  the  solution  bl9^a  \q).  Bie  first 
eritleal  point  encountered  is  the  required  critical  length. 

In  the  second  approach,  we  solve  a  systea  of  linear  differential 
enuatlons  with  a  tvo-polnt  boundary  condition,  “ails  aethod  is  closely 
allied  to  classical  eigenvalue  techniques,  with  the  length  x  playing 
the  role  of  an  eigenvalue. 

In  the  single  ti»e- Independent,  energy-independent  case  treated 
above,  both  ■ethods  are  readily  pursued,  both  yielding  explicit  analytic 
expressions.  It  is  only  when  we  turn  to  the  discussion  of  nore  realistic 
proossses  that  a  fork  In  the  road  develops.  A  foimdatlon  which  yields 
linear  ecuatlons  appears  very  nvieh  nore  desirable  If  a  theoretioal 
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traatanit  is  oontfl^pl&t«d,  vtilla  an  approach  buod  v^oc  nonltnaar 
eqaatlOM  dvtaraljwd  b/  laltlal  oondltloDS  appoara  prafarabl*  if  a 
maarleal  traat— nt  is  dsslrsd. 

In  futurs  vork,  vs  shall  discuss  tbess  points  In  grsatsr  detail. 

6.  Intarllnking  of  tbs  Two  Tssfanlquss. 

BsTlncs  seen  that  tvo  qnlts  different  analytic  approaobss  to  tbs 
studjr  at  navtron  transport  are  possible,  vs  visb  to  detenlne  tbe 
rslatiOBS,  if  any,  that  exist  aaoDg  tbe  various  fimetioDS  vs  ba-vv 
defined.  Ihls  Unirtafl  of  teebniques  can  be  useful  in  tvo  ways.  On 
tbe  oca  hand,  it  sbova  us  hov  to  linearise  oertaln  types  of  Rioeatlan 
functional  equations,  and  on  the  other  band,  it  sbovs  hov  to  replace 
two-point  bovBidary  value  problans  involving  linear  equations  by 
initial  value  problens  for  nonlinear  equations. 

In  order  to  obtain  the  desired  rslatiens,  ve  cosine  the  tvo 
fonmlations  in  tbe  foUwlng  Banner.  In  calculating  qp(n,y),  ve 
first  observe  that  as  a  result  of  one  neutron  entering  [7,0]  froa  tbe 
left,  vhieb  is  to  say  leaving  [^x,y  j  frea  tbe  right,  ve  have  a 
probability  p(k,y),  in  tbe  notation  of  Section  S,  of  having  t  neutrons 
reflected.  lisre  we  are  speaking  only  of  those  neutrons  that  eaerge 
froa  [y,0  j  ,  aovlng  to  the  left,  for  tbe  first  tiaa.  Sach  of  these  k 
neutrons  can  now  be  coneidered  as  the  initiator  of  a  new  process 
starting  froa  y  in  which  a  nsutron  enters  [x,y3  tram  the  right,  tbe 
particle  aovlag  to  the  left. 

In  these  texas,  let  us  see  bow  ve  can  obtain  exactly  n  neutrons 
eaarging  froa  x  over  all  tlas.  Suppose  that  the  first  of  the  neutrons 
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daserlbed  ab<rr»  produces  osutrone,  the  secood  n^,  sod  the  k-th,  a 
quaxitlty  vbsre 

(6.1)  ♦  Og  ■*•...  -  n. 

Since  the  probability  of  this  event  Is 

(6.2) 

ve  derive  the  reliitlon 

(6.5)  qjj(n,y)  =  £  p(k,y)  f  Z'ljoi.y)  ' 

”  k«l  *■ 

where  the  inner  svMatloD  is  taken  over  the  lsttle»-points  satisfying 
eqxntlon  (6.1). 

ReosLlllng  the  definitions 

52  *•“  -  g(r,y), 

n-0  ” 

(6.U)  Z  -  f(r,y), 

DoO 

52  p(o,x)  r“  o  h(x,r), 
n^O 

ve  oi>taln  froB  equation  (6.5)  a  relation  cooneoting  the  generating 
functions, 

(6.5)  g(r,y)  =«  h(x,f(r,y)). 

An  aanlogous  relation  ean  he  derived  if  ve  consider  q^(y,r) 


Initially. 
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tteM  r«l«tlaDS  ve  mux  terlv*  •quatlOD*  eaDnaatine  tba  axpaetad 
ral«**.  Howarar,  as  aa  ahall  aaa  balow,  tbasa  aquation*  oan  readily  b* 
obtalxiad  dlraetly. 

•Om  I  — liiilT  of  the  papar  all!  be  devotad  solaly  to  a  (U*aua*lon 
of  axpaatad  Talae*,  Blnaa  tha  aaalacua*  of  the  result*  abora  for 
gaaeratias  fiaetloa*  are  man  ralatlon*  eanaaetlng  aharaetarlstia 
fiawtlonal*.  aloe*  eonsldaratloa  of  tha**  anti  tie*  aould  both  et»- 
pUeat*  th*  saalyai*  and  alevat*  tha  ocnoaptual  level,  w*  have  kept 
thaa  offsta^  until  aaothar  data. 

7.  Kqeetad  Value*. 

Let  vm  now  ladloata  how  to  obtain  th*  aquation*  for  axpaetad  value* 
vlthocct  throacb  tha  Intezaedlats  step  of  prohahillty  dlstrlbutloa* . 

Ihl*  dlseuaslon  is  preparatory  to  the  oour**  we  shall  pursue  In  th* 
traataenk  of  the  aora  eoivlleatad  tlae-dapandent  and  enargy-dapendant 
eases.  V*  any  of  tha  proeasse*  a*  being  dataxalnlstla,  a*  Is 

coHHB  In  radlatlw*  transfer  theory. 

Let  U(x)  danota  th*  axpaatad  m^ber  of  neutrons  refleetad  trm 
[  x,o3  over  tlae  a*  the  result  of  an  inddant  trigger  neutron  at 
tlae  laro.  In  passing  through  th*  sa^nt  fx,*  -  A]  th*  trigger 
neutron  ney  mdargo  fission.  In  which  ensa  one  neutron  Is  raflaetad 
and  one  1*  Incident  on  [x  -  A,  0]  ,  the  expected  n\ad)*r  of  reflacted 
neutrons  frca  this  secant  being  U(x  -  A).  If  the  trigger  neutron 
doe*  not  undergo  fleslon,  then  th*  axpaetad  nvahar  of  neutrons 
rafleatad  dlraetly  1*  U(x  -  A).  In  addition,  one  of  tha  neutron* 
reflacted  fron  the  lepwnt  [x  -  A,  o]  nay  undergo  flsalon  giving  rise 
to  another  trigger  nsutron  Ineldant  on  ths  sepMnt  [x  -  A,  0]  ,  which 
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In  turn  gives  rise  to  the  expected  nvaber  of  reflected  neuttoxw 
U(x  -  d).  All  other  prooesses  lead  to  terse  of  at  least  second  order 
in  d.  Vedcing  account  of  these  possibilities,  ve  obtain  the  equation 

(7.1)  U(x)  -  x‘^d  (1  +  U(x  -  d))  +  (1  -  A'^d)  U(x  -  d) 

♦  (1  -  d)  U(x  -  d)  X‘^  dU(x  -  d)  ♦  o(d). 

Letting  d-e  0,  ve  obtain  the  differential  equation 

(7.2)  U'(x)  -  [l  +  U^(x)]  ,  U(0)  -  0. 

In  addition,  if  ve  let 

V(x)  -  the  expected  nimber  of  neutrons  eaergcnt  fr«  one  end 
of  a  rod  of  length  x  as  a  result  of  a  trigger  neutron 
Incident  at  the  other  end, 

ve  find  that  v(x)  satisfies  the  equation 

(7.3)  V(x)  =  U(x)  V(x),  V(0)  .  1. 

Smilarly,  ve  see  that  «p(y)  and  Uj^(y)  satisfy  the  linear 
e  ^uatlons 

u^(y)  =  x"^ 

(7.U) 

u£(y)  =  -X"^ 

To  derive  equations  (7.1^)  ve  trace  the  trigger  neutron  through 
either  the  se^est  [y  ♦  d,y]  or  [y,y  -  d  ]  ,  as  the  case  aay  be. 
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Coasldftrlac  thu  oam  la  vhlcb  It  aovas  to  the  rlgbt,  w«  find  that  If  It 
doas  not  uadarso  flaaloa  on  ttaa  aof^at  [  7  ♦  A,y  ]  ,  la  affect  va  have 
a  proeaaa  la  vhldi  there  la  a  trigger  nautroc  aovloc  to  the  right  at 
y.  If  It  doea  tBdarco  flaaioa  than  va  have  a  prooaaa  In  vhleh  there 
are  trlxgar  aeatroBa  Borlag  tovard  tbs  right  and  toward  the  left  at  y. 
Mtlag  Into  aeooiat  the  frohahlUtlea  of  these  e'/enta  va  have 

(7.5)  ttpCy  ♦  A)  -  (1  -  ^*^A)  Up(y)  ♦  X"^A^Up(y)  +  u^(y)  j*  o(a). 
Slailarly  for  u^(y)/  ve  find 

(7.6)  Uj^(y)  -  (1  -X"^A)  Uj^(y  ♦&)■*.  X’^A  ♦  A)  ♦  1^(7  +  A)|  +  o(a). 

BituBtioBs  (7.^)  follow  directly  froB  these. 

FoUowlag  the  argvBsnt  given  In  the  previous  section,  it  is  easy 
to  see  that  ve  obtain  the  foUcadag  two  equations 

V»>  -  W  \ij), 

(7.7) 

«j^(y)  -  V(x  -  y)  ♦  U(x  -  y)  ^{7), 

vhldi  eaoaeet  the  solutions  to  the  linear  and  nm^lnesr  differential 
equations  given  above. 

Consequantly  the  tranafomatlon  fozmilas  are 

' 

(7.8) 


and 
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(7.9) 

V(x.y)  3  Uj^(y)  -  uj,(y)  . 

Mans  of  those  relatlOBs,  «e  can  pass  froi  tbs  solution  of 
one  type  of  equatlcn  to  the  other.  We  have  thus  obtained  a  set  of 
rslstloos  linking  the  two  types  of  functions. 

S.  Tlne-and-Kaergy-Independeat,  Variable  Observer. 

Let 

Vj^(z;x)  -  the  expected  niaber  of  neutrons  passing  z  toward 

the  right  over  all  tljse  as  a  result  of  one  trigger 
neirtron  cnterlUt;  'the  bar  [x,0]  at  the  left,  at  x, 

v^(z:x)  >  the  corresponding  nvaiber  passing  to  the  left. 

Then  ve  easily  derive  the  eqiiatlons 

Wj^(z  +  A;x)  =  Vp(m)x)  -  Wj^(z;x)  -►  o(a), 

(3.1) 

Wj^(r  +  A;x)  -  Vj^(z;x)  +  Vj^(z;x)  +  o(a), 
which  lead  to  the  linear  systen 
w^(z;x)  =  -  X“^wj^(z;x) 

(8.2) 

w£(z;x)  =  Wj^(z;x), 

where  the  prlass  denote  differentiation  with  respect  to  z.  The  two 
point  boundary  eoadltioos  are 
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(8.3)  Wp(x;x)  -  1,  Wj^(Oix)  .  0. 


Co^urlag  tbcM  ei^uatlan*  and  bousdary  eonditlona  vltb  those 
obtained  for  Up  and  la  the  previous  section,  ve  find  that 


Uj^(l)  -  Vjj(tix), 

(S.U) 

Ujj(t)  -  Wj^(t;x) . 


This  systea  constitutes  a  reciprocity  lav  vhleh  states,  In  part,  that 
the  expected  maiber  of  neutrons  eaerglng  fron  x  over  all  tine  as  the 
result  of  a  trigger  neutron  aovlng  to  tbs  left  at  x  Is  equal  to  the 
expected  mssber  of  neutrons,  over  all  tine,  aovlng  to  the  right  at  x 
as  a  result  of  a  trigger  neutron  Introduoed  at  x.  'Ilw  reciprocity 
obviously  bolds  under  nueb  nore  general  elrcuutances . 


9-  The  Internal  flux . 

Let  us  nov  consider  the  expected  mssber  of  neutrons  aovlng  to 
the  right  which  pass  a  fixed  point  x  as  a  result  of  a  source  neutron 
aovlng  to  the  right  at  a  fixed  point  y,  y>x  (see  Pig.  3).  Ve  consider 
this  as  a  function  of  the  length  of  the  bar,  x,  say  f(x:y,x),  y  and  x 
being  held  fixed.  This  will  Illustrate 


Source  Obsez^r 


Figure  3 


how  Internal  flows  aay  be  deteralnad.  If  x  ay,  then 
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(9.1)  r(x;y,z)  -  f(yiy,t)  -  Vj^(i;y) 

In  the  notatlQB  of  th*  prsvloiu  Mctlon.  If  x>y, 

(9.2)  f(x;y,r)  -  Vj^(x;y)  ♦  tj(y)  u(x-y)  Vj^(t;y). 

■Rje  flrrt  tem  in  the  right  hand  side  of  the  above  eqtMtlon  la  the 
direct  flux  and  the  second  la  the  secondary  flux  dite  to  reflections 
fro*  the  segBents  [y,©]  aad[x,y  j  . 
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III.  Tla»-Depm><ttnt ,  Saergy-lndepmdeBt  Ca— 

10.  Expected  Valu». 

Let  us  Dov  oonslder  a  tlae-depandent  prooass.  It  was  pointed  out 
to  the  authors  by  T.  B.  Barrls  that  the  teehniquss  used  in  the  preceding 
seetioBS  oan  be  equally  well  e^plosred  in  the  diseuaslon  of  tljee* 
depeadenee,  provided  that  wa  Introduce  the  iaportant  function 
prtAiabllity  density. 

We  define  the  follovlng  two  fisietlona. 


(10.1) 


U(x,t)  «  the  expected  n»ber  of  particles  easrglng  fros  x 

between  tins  sero  and  tine  t  frcB  a  rod  of  length  x 
due  to  a  trigger  neutron  entering  x  at  tiai  sero. 

u(x,t)dt«the  expected  niaber  between  t  and  t  4-  dt  arising 
In  the  easts  Banner. 


Rote  that  •  u. 

To  obtain  an  equation  for  U(x,t),  we  begin  with  the  fact  that  a 
neutron  traversing  [x  4-  d,x]  splits  with  probability  x~^A  and  continues 
unaltered  with  probability  (l  -  A  ^h).  Ibe  neutron  which  has  not 
suffered  a  collision  produces  a  certain  flux  of  neutroxis  fron  [x.ol 
In  the  tljns  Interval  [x,s  ♦  ds]  ,  naasly  u(x,a)  ds.  'Rie  velocity  of 
a  neutron  Is  denoted  by  o. 

Of  this  flux,  an  expected  aaount  (l  -X'^A)  u(x,s)  ds  continues 
through  *  d,x3  without  eoUlslon,  while  A*^Au(x,s)  ds  produces 
neutrons  going  In  both  dlreotlons.  Ihe  quantity  a~^Au(x,s)  ds  going 
to  the  right  produces  an  expected  reflection  of  U(x,t  -  s)  ‘>A~^Au(x,s)  ds, 
during  the  reaalnii^  tlas. 
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Addlnc  the8«  effect*  together,  ve  see  that 

U(x  +  h,t  1^)  -  [U(x,t)-.l]  4.  (1  -  X’^A)  |^U(x,t) 

■*•  f  u(x,s)  X'^A  U(x,t-s)  ds  I  *  o(a). 

0  i 

Passing  to  the  liMlt,  -e  obtain  the  Integro-differentlal  equation 

t 

(10.3)  ^x  c  ^t  ’  ^  ^  “(x#*)  U(x,t-8)  ds  +  x”^. 

0 

'n^e  initial  conditions  in  space  and  tlae  are 
U(x,0)  ^  0, 

(10.1*) 

U(0,t)  -  0. 

Using  Laplace  tranafom  techniques,  we  can  solve  for  U  since  the 
equation  in  (10.3)  is  of  convolwctlon  type.  A  nvad>er  of  interesting 
quest  ions  arise  in  this  analysis  which  we  shall  discuss  at  a  later 
date. 

Let 

V(x,t)  «  tile  expected  nmber  of  neutrons  transaltted  through 
a  bar  of  length  x  in  tiae  t. 

Then  we  obtain  the  equation 

t 

(10.5)  V(x*.A,t  f )  -  V(x,t)  4.  (1  -  X'^A)  ^  u(x,*)  x'^A  V(x,t-*)  da  4-  o(A) 

0 
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which  leads  to  the  Intagro-dlffcrcntlal  equation 

t 

(10.6)  '^x  *  \  “(*»•)  V(x,t-a)  da. 

0 

In  addition  V  aatlsfies  the  conditions 
V(0,t)  -  1 

(10.7) 

V(x,0)  •  0. 

11 ■  Hovlne  Source . 

Let  us  now  describe  the  process  In  terns  of  a  variable  source  let 

u^(yjt)  ■  the  expected  ntaber  of  neutrons  energlng  fr<ni  x 
over  the  tli»  interval  [0,t]  due  to  an  initial 
neutron  novlng  to  the  left  at  y  at  tine  sero. 

Up(y,t)  -  the  corresponding  quantity  due  to  a  neutron  ■O'.'lng 
to  the  right. 

then  we  readily  obtain  the  following  equations: 

Uj^(y,t)  »  *  ti,  t  -  ^)  *  Ujj(y  ♦  A,  t  -  -)^ 

♦  (1  -  x”^A)  Uj^(y  A,  t  -  ^)  ♦  o(a). 

(11.2) 

Uj,(y,t)  -  x'^A  |’u^(y  -  A,  t  -  |)  ♦  Uj^Cy  -  A,  t  -  f)^ 

♦  (1  -  x"^A)  Uj^(y  -  A),  t  -  |)  ••■  o(A). 
rren  whence  we  obtain  the  systsn  of  differential  equations 
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(11.5) 


1  S  S 


X 


-1 


\ 


for  C<y<x,  with  the  coaditiocis 


Uj^(y,0)  -  0,  \(*>^)  -  1. 

(11.^) 

Up(y,o)  -  0,  Up(o,t)  -  0. 

Proceeding  as  in  the  stationary  case  we  readily  obtain  the  follcrwlng 
relations,  shot. log  the  Intercodnectlons: 

t 

Vy,t)  -  \  u(y,8)  Uj^(y,t-8)  ds, 

0 

(11-5)  t 

Uj^(y»t)  *  V(x-y,t)  4.  ^  u(x-y,B)  Up(y,t-8)  ds. 

0 


IV.  Tl«e  lDdapM>d«nt-teTgy  Dependent  0>— 


12 .  Ixpeeted  Valuas . 

Ltt  us  BOV  ecxisldtr  the  eas«  wttare  va  take  Into  aceount  ths  fact 
that  the  Teloeltj,  or  emrgy,  of  tbs  partlele  detsralasB  ths  probability 
of  fission.  FurtherBora,  fission  vlU  produce  tvo  partlelas  with 
energies  different,  in  general,  froa  that  of  the  parent.  For  con* 
venlance,  0  and  x  in  Fig.  1  are  interchanged: 


>■  > 

0  d 


Flgxire  4 


We  define  the  functions 


u(x,v;r)  dr  «  the  expected  nuaber  of  neutrons  reflected  frcB 
C  with  energy  betveen  r  fr»d  r  ♦  dr,  over  all 
tlae,  due  to  an  Initial  neutron  of  energy  v 
entering  the  rod  at  0  at  tlae  aero. 

(12.1) 

v(x,v;r)  dr  >  the  expected  nvaiber  transaltted. 

4(s,r)  dr  •  the  probability  of  energy  change  froa  level  s 
to  a  level  betveen  r  and  r  -f  dr  due  to  a 
fissioning. 

®(^)  ■  probability  of  a  fissioning  in  a  segpent 
of  length  d  of  a  neutron  at  energy  level  v. 
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u(*,wjr)  -  (1  -  7-^)  u(x  -  A,wjr)  (1  -  — ) 
(12.2)  ♦  (1  -  x(v)  )  3  ^  u(x  -  q(i,r)  -yrj^ 


(1  -  xfv)  ^  ^  ^  ^  “(*  "  3,v;b)  u{*  -  A,t;r) 


0  0 


-5^-^  q(v,v)  1.  ^  u(*,g;r)  q(w,i)  d*l 

*■0  J 


PuslAg  to  the  Halt,  ve  obtain  the  integro-dlfferentlal  equation 


(12.3) 


du 

’S 


da 

Aliy 


+  )  )  q(»,t)  u(x,w;a)  u(x,t;r) 
0  0 


*  xiry  ^  “(*»*;>•)  q(v,*)  d*. 

0 

Ikera  la  a  alallar  equetloa  for  V,  the  tranaalttod  flux. 

Let  ua  now  obtain  the  funetlona  whloh  will  peralt  xia  to  llnaeLrlte 
theae  equatloDa.  Introdiwa  tha  funetlona 
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>l|(7 ^  c(p«et*d  Biaiber  of  neutroDa  — rglng  froa  0 

vitb  anargsr  batiwn  r  aad  r  •«■  4r  ovar  all  tlaa 
Ana  to  aa  laitlal  DautroD  of  anargy  v  atartiag 
(12.^)  to  ttaa  right  at  j,  0<y<x,  at  tlaa  aaro. 

^  ~  eorraapondlng  quantity  vttan  tha  nautron 
atarta  to  tha  laft. 


f  aa  bafora, 


»L(7,»ir)  -a- 


^  Up(y,a;r)  q(a,a) 


(12.5) 


•-  0 


♦  ^  Uj^(7»»;r)  q(*,»)  4a 


♦  o(a). 


vith  a  sorraapcDdlng  aquation  for  Up(y,v)r). 

To  ohtala  aquatlona  eonaaatlng  thaaa  fmetlona,  wa  again  prooaad 
aa  In  tha  pravloua  aaetlona.  tha  raault  la 


«j,(7»v}r)  -  ^  u(x  -  y,w;a)  •  Uj^(y,a;r)  Aa, 

0 

-  ▼(y.wjv)  ♦  ^  u(y,w;a)  Uj^(y,a;r)  Aa. 


(12.6) 
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